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ABSTRACT 

We re-examine the evidence for the existence of ultra-massive (M > 1.1M Q ) white 
dwarfs based on gravitational redshift of white dwarfs in common proper motion bi- 
naries or in clusters, on parallax measurements, on orbital solutions, and, finally, on 
the analysis of hydrogen line profiles. We conclude that the best evidence is largely 
based on the analysis of Balmer line profiles although the companion to the A8V star 
HR 8210 is a compelling case made initially using the large binary mass function and 
confirmed by an analysis of the Lyman line spectrum. The confirmation and identifica- 
tion of high-mass white dwarfs, more particularly non-DA white dwarfs, using parallax 
measurements may prove critical in establishing the population fraction of these ob- 
jects and in constraining the high-end of empirical initial-mass to final-mass relations. 
The existence of a substantial population of ultra-massive white dwarfs supports the 
concept of a steeper initial-mass to final-mass relations linking 6 Mq progenitors with 
> 1.1 Mq white dwarfs. 
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• ' 1 INTRODUCTION 

Ultra-massive white dwarfs, generally denned as having 
masses exceeding 1.1 Mq, remained relatively rare un- 
til deep spectroscopic surveys uncovered these intrinsically 
faint o bjects. The hot hyd rogen-rich (DA) white dwarf 
GD 50 dBergeron et al l 1 199 j) . the carbon-rich white dwarf 
G35-26 (|Theill et all Il990t ), and two DA white dwarfs 
(PG 165 8+441 and PG 0136+ 251) from the Palomar-Green 
survey (|Schmidt et all 1 19921) were rare examples of this 
phenomenon. Interestingly, spectroscopic follow-up of ex- 
treme ultraviolet (EUV) surveys of the local (d < 100 pc), 
hot white dwarf population managed to identify many new 
ultra-massive white dwarfs on the basis of a large, spectro- 
scopically determined surface gravity (loggr > 9). Based on 
the Extreme Ultraviolet Explor er (EUVE) and ROSAT Wid e 
Field Camera (WFQ surveys, IVennes et all (Il996l. Il997bl) 
IVennesI i| 19991) . iFinlev et all <| 19971 ). and lMarsh et all (|l997l ) 
added a dozen new objects to the population. More recently, 
a re-an alysis of the Palomar -Green (PG) sample of DA white 
dwarfs ijLiebert et alll2005l ) and a study of the white dwarf 
mass distribution in the Sloan Digital Sky Survey (SDSS; 
Kepi er et al . 2007) seem to generate similar yields of ultra- 
massive white dwarfs. 

Are these objects products of single star evolution or 
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double degenerate mergers? IWeidemannl (|200d ) reviews the- 
oretical arguments in favour of an upper mass limit for 
white dwarf stars larger than the canonical upper limit of 
1.1 Mq, and possibly as large as ~ 1.3 Mq. A limit of 
1.1 Mq is generally assumed because carbon ignition in high 
mass cores would lead to thermonuclear runaway or core 
collapse. However, the effect of mass loss may alter the sce- 
nario and lead to the fo rmation of a massive oxygen-neon- 
magne sium white dwarf ( Nomotolll984l ) . iGarcia-Berro et all 
ljl997h and iRitossa et al. ( 19961 ) successfully evolved 9 and 
10 Mq stars, which, following an off-centre carbon ignition 
in partial electron degenerate conditions, generated oxygen- 
neon core white dwarfs with carbon-oxygen shells, and to- 
tal masses of 1.15 and 1.26 Mq, respectively. On the other 
hand, the mer ger scenario proposed f or the origin of Type 
la supernovae (Yungclso n et al] 1 19941 ) also generates ultra- 
massive white dwarfs. However, the merg er process itself 
and t he fate of these objects are uncertain (|Segretain et all 
I1997D . 

We re-examine the evidence for the existence of ultra- 
massive white dwarfs. Mass measurements are based on 
gravitational redshift measurements (§2), radius (paral- 
lax) measurements (§3), orbital parameters (§4), and sur- 
face gravity measurements (§5). The gravitational redshift 
measurements (tx M/R), radius measurements, and sur- 
face gravity measurements (oc M/R 2 ) are converted into 
mass measurements by adopting mass-radius relations for 
a variety of model structures. We adopted the models of 
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Table 1. The CPMB Sample" 



WD CPMB Types 



0433- 


h470 


LP 356-525/BD+26 730 


DA +dK5 


0738- 


-172 


LP 783-3/LP783-2 


DZAQ+dM6 


0820- 


-585 


LP 186-119/LP186-120 


DA+dM3 


1750- 


^098 


G140-B1B/G140-B1A 


DC+dK2 


2323- 


-241 6 


G275-B16A/G275-B16B 


F7+dK7 



" ISilvestri et al.l l|200ll l. 

b Not a white dwarf, see text. 



iBenvenuto fc Althausl ij 19991 ) with a C/O core and zero 
metallicity, a helium mantle (logg(He) = —2), and with a 
thick (log g(H) = —4) hydrogen layer to describe hydrogen- 
rich DA white dwarfs or without a hydrogen layer to describe 
non-DA white dwarfs. In §5, we also present evidence that 
deep surveys of high-proper motion white dwarfs are now 
extending the ultra-massive population toward lower tem- 
peratures, and it is estimated that w 10% of white dwarfs 
have masses in excess of 1.1 Mq. We summarize and con- 
clude in §6. 



2 BASED ON GRAVITATIONAL REDSHIFT: 
THE CPMB SAMPLE AND CLUSTER 
MEMBERS 

The mass ive white dwarf an d member of the Pleiades 
LB 1497 (We gner et al.l Il99ll ) remained a unique candi- 
date for having a mass in excess of 1.0 Mq and possibly 
in excess of 1.1 Mq among stars with gravitational redshift 
measurem ents dWegner fc Reidl Il99ll : iBergeron et al.l 1 19951 ; 
lReidlll996l 'l . However, applying the same tech niques to white 
dwarf s in common proper motion binaries ISilvestri et al.l 
|200ll ) proposed new ultra-massive white dwarf candidates 
(Table [TJ. We demonstrate that this is not the case and 
that the ultra - massiv e white dwarf sequence presented by 
ISilvestri et all (|200lT ) and based on gravitational redshift 
measurements does not exist. We also re-evaluate the mass 
of LB 1497 and show that its case remains compelling. 

2.1 WD 0433+470: LP 356-525 

The hydrogen-rich (DA) white dwarf WD 0433+470 is part 
of a CPMB with the putative Hyades member BD+26 730. 
However, the presence of a cool, presumably old white 
dwarf in the Hyade s poses a chronology problem that 
ICatalan et alj l|2008bh tentatively resolved by substituting a 
normal C/O core with a Fe core thereby shortening the cool- 
ing age of the white dwarf from 4 Gyr to lG yr. Their analy- 
sis is b ase d on the Balmer line pro file fits of Ber geron et al.l 
|200ll ) and lZuckerman et alj (|2003l ) . Adopting a straight av- 
erage of these measurements, T c g = 5527+130 K and loggr = 
8.08 + 0.10, we determined a mass M = 0.63 + 0.06 Mq and 
absolute visual magnitude My = 14.78 ± 0.30 us i ng th e 
DA mass-radius relations of IBenvenuto fc Althausl (l999). 
The resulting mass is well below th e mass based on the 
gravitational redshift measurements of lSilvestri et al.1 (|200ll ) 
M = 1.12 ± 0.01 Mq. Moreover, the apparent (V = 15.81) 
magnitude combined with the absolute visual magnitude for 
the low-mass solution corresponds to a distance modulus of 



m - M = 1.03 ± 0.30 (d = 16 ± 2 pc) in agreement with 
the Hipparcos parallax measurement for BD+26 730, it = 
56.0 ± 1.2 mas (jGould fc Chaname|[200i ) or d — 17.9 ± 0.4 
pc, a nd the white dwarf parall ax measurement, ir = 60 ± 3 
mas (lHarrington fc Dahnll980h or d = 16.6+0.8 pc. Accord- 
ing to lCatalan et al.l (|2008bh Tthe mass of an iron core white 
dwarf would be lower by ~ 0.2 Mq for the same temper- 
ature and surface gravity. The white dwarf radius inferred 
for this configuration would also decrease by 18% resulting 
in a closer distance for the white dwarf d = 13 ± 2 pc that 
is marginally consistent with the parallax measurements. 

2.2 WD 0738-172: LP 783-3 

The metal-rich hydrogen-poor white dwarf LP 783-3 
(DZAQ6) is paired with the dM sta r LP 783-2. The paralla x 
of the white dwarf, it — 112+3 mas l|van Altena et al]ll995l ). 
translates into a distance modulus of m — M — —0.25 ± 0.06 
and, using V = 13.03, we determined the absolute magni- 
tude M v = 13.28 ± 0.07. 

By fitting optical/infrared spectral energy distribution 
and by co nstraining the absolu te magnitude to Mv = 
13.31 ±0.05 lBergeron et""aH (|200lh determined T cB = 7710 + 
220 K and log g = 8.09 ± 0.03. 

We revisited this result by computing a H/He model 
grid for effective temperatures 7000 ^ T c g ^ 8500 K 
and surface gravities 7.0 ^ log<? ^ 8.5, and a hydro- 
gen abundan ce log H/He = —3 co nsistent with the mea- 
surements of IBergeron et all (|200ll ). The absolute magni- 
tudes were calculated using the mass radius relations of 
IBenvenuto fc Althausl (1 19991) for non-DA stars. We adopted 
the temperature of lBergeron et all <|200lT) . T eS = 7710+220, 
and the absolute magnitude Mv = 13.28 ± 0.07 and prop- 
agated the errors. We determined a surface gravity loggr = 
8.12+0.09 that corresponds to a mass ofM = 0.64+0.06M©. 
The mass of the white dwarf is normal. 

Adopting the distance modulus for the white dwarf 
(m — M = —0.25 ± 0.06) we estimate an absolute lumi- 
nosity for the red dwarf companio n My = 16.92 ± 0.06 that 
corre sponds to a spectral dM6.5 (|Kirkpatrick fc McCarthy! 
Il994h . 

2.3 WD 0820-585: LP 186-119 

The DA white dwarf (L186-119) is paired with a dM star 
(L186-120). We obtained a low- dispersion (w6 A) EFOSC2 
spectrum from the ESO archive^] and fitted the Balmer line 
profiles from H/3 to H9 u sing our most recent pur e hydrogen 
model atmosphere grid (|Kawka fc Vennea 120061 ') . Figure Q] 
shows the results of our analysis. Using the DA mass-radius 
relations of lBenvenuto fc Althausl I\199$ ) , the surface gravity 
and effective temperature measurements (log g — 8.12+0.08, 
T cft = 9310 ± 70 K) correspond to a mass M = 0.67 ± 
0.05 Mq well below the mass of M = 1.10 ± 0.03 M Q esti- 
mated by ISilvestri et all (1200 ll ) and based on their gravita- 
tional redshift measurement. The calculated absolute mag- 
nitude is My — 12.66 ± 0.22 and corresponds to a distance 
modulus m - M = 3.9 ± 0.2 (d = 60 ± 6 pc). 

1 Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under programme 078.D-0824. 
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Figure 1. Model atmosphere analysis of the DA white dwarf 
WD0820-585 (L186-119). 



ISilvestri et ail (|200lh listed the spectral type of the 
companion as dM5. We revised this estimate using their 
V photometry as well as 2MASS JHK measurements 
i|Cutri et alj l2006h . The measured V - J = 4.0 index 
corresponds to a »dM3.5 (M y = 12.05) spectral type 
|Kirkpatrick fc McCart"hvlll994l ) and the corresponding dis- 
tance modulus (m — M = 3.8) confirms the physical associ- 
ation of the two stars and the normal mass (M « 0.7 Mq ) 
for the white dwarf component. 



2.4 WD 1750+098: G140-B1B 

The DC white dwarf G140-B1B is paired with the K star 
G140-B1A (HD 162867). The optical/2MASS colour index 
V — J = 1. 68 ± 0.08 of G140-B1 A constrains the spectral 
type to K2 (|Bessell fc Brettl 19881) and the absolute magni- 
tude to My - = 6.2 ± 0.2 | |Dotter et al.ll2008l ). Therefore, the 
distance modulus is m — M — 3.2 ± 0.2 and corresponds to 
a dista nce of d = 44 pc. Adoptin g V = 15.72 for the white 
dwarf l|Eggen fc Greensteinl Il96a ) . the absolute luminosity 
is My = 12.5 ± 0.2. To determine the corresponding stel- 
lar parameters we computed a grid of pure He models at 
9000 < T cff < 11000 K and surface gravities 8.0 < \ogg < 
9.0 and calculated the absol ute magnitudes using the non - 
DA mass-radius relations of iBenvenuto fc Althaud (j 1998T) . 
Combined with the colour index V - I = 0.178 ± 0.060 
ijSilvestri et alj|200ll ). the absolute magnitude corresponds 
to T eff = 9200 ± 600 K and logg = 8.14 ± 0.25 that trans- 
late into a mass of M = 0.66 ± 0.16 Mq. T he presence of 
Ha in the white dwarf spectra is uncertain as lSilvestri et al.l 
J200ll) andlWegner fc Reidl (|l99ll ) claim Ha detections while 
iReidl (fl996l ) does not report the detection. The absence of 
Ha would invali date the white dwarf gravitational redshift 
measurement of ISilvestri et all (|200ll ). However, the coin- 
cident proper motions and the normal mass for the white 
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Figure 2. ESO UVES spectrum of the F7 star G275-B16A 
(shifted up by 0.5 in relative flux units) compared to low- 
metallicity models, [M/H]= -1.5 and [M/H]= -2.0 (shifted up 
by 1.0 in relative flux units). 



dwarf, which implies similar distance moduli for the two 
stars, strongly favour a physical association for the pair. 



2.5 2323 241: G275-B16A 



ISilvestri et alj (|200ll ) estimated a white dwarf gravitational 
redshift v g = 132 ± 1 km s^ 1 that translated into a mass 
M — 1.19 ± 0.01 Mq. This conclusion is based on the as- 
sumption that the star is a white dwarf part of a common 
proper motion binary and that the velocity differential is en- 
tirely caused by the gravitational redshift of the white dwarf. 
The companion is a late type star. We now re-examine the 
evidence. 

First, Figure [2] shows a comparison between two model 
spectra and the UVES spectrum from the ESO archive^ 
of G275- B16A. The model sp ectra were obtained from the 
library of iMunari et alj (|2005r Fl which is based on Kurucz 
models. Models and data are presented at a resolution of 1 A. 
The Balmer line series clearly extends to H15 and implies 
a normal main-sequence surface gravity and a low metallic- 
ity. Indeed, we estimate [M/H]~ -1.5 and T cS = 6500 K 
(log g = 4.5). Therefore, the star G275-B16A is not a white 
dwarf but a low metallicity main sequence star. 

Next, the distant modulus and proper motion measure- 
ments suggest that the pair is coinci dental. Based on 2M ASS 
measurements and V magnitudes (jSilvestri et alj [20051 ) we 
estimate V — J = 0.99 and 2.43 for G275-B16A and G275- 
B16B, respectively. The colou rs correspond to spectr al types 
of F7 and K7, respectively l|Bessell fc Brettl Il988l). Using 
the low metallicity ([M/H]« —1.5) models of iDotter et alj 

2 Based on observations made with ESO telescopes at the La 
Silla Paranal Observatory under programme 165.H-0588. 

3 Available at |http://archives. pd. astro. it/2500-10500/| 
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(2008) for the F star, the colour index V - J = 0.99 corre- 
sponds to an absolute magnitude My = 5.2 (or T e « = 6200 
K). A somewhat brighter absolute magnitude My = 4.9 
is obtained with the model at T e fi = 6500 K. Note that 
low-metallicity stars are fainter than their solar-metallicity 
counterparts. Using the solar-metallicity models for the K 
star we determined a temperature of T c s = 4200 K and an 
absolute magnitude My = 8.0. Therefore, the F star is at a 
considerably larger distance (m — M — 10.6-10.9, d — 1.3- 
1.5 kpc) than the K star (m — M = 7.7, d = 350 pc) and 
the pair is purely coincidental. 

This conclusion is supported by proper motion and ra- 
dial velocity measurements tha t place the F star on a dif- 
ferent course than the K star. IZacharias et al.l l|2005l ) list 
the proper motions (in mas yr^ 1 ) as /i Q = ± 3 and 
/i<5 = —36 ± 5 for the F star, fi a = —5 ± 5 and us = —64 ± 5 
for the K star. Consequently, the di fference in velocities , 
w(F star) - w(K star) « 132 km s" 1 l|Silvestri et al.ll200ll ). 
may no longer be interpreted as being due to the effect of 
gravitational redshift, but as being due to diverging space 
motions. 

The F star space motion is c haracteristic of low- 
metallicity stars (|Chiba fc Beerslf2000t ) . Adopting a distance 
of 1.4 kpc the radial velocity (f ra d = 133 km s -1 ) and 
prope r-motion correspond, following Ijohnson fc Soderbloml 
1 19871 ). to a Galactic velocity (in km s" 1 ) vector (U, V, W) = 
(143, —185, —138), which confirms its membership into the 
old disk or possibly the halo. 



2.6 WD 0349+247: LB 1497 

LB 14 97 is a member of the Pleiades and IWegner et alj 
|l99ll ) reported a gravitational re dshift of v„ = 8 4 ± 9 
km s -1 for this DA white dwarf. iBergeron et al.l (|l995l ) 
also resported a temperature T e g = 31660 ± 350 K and a 
high gravity of log o = 8.78 ± 0.05. Usin g the DA mass- 
radius relations of Benvcnuto & Althaus (1999) and the 
temperature from IBergeron et alj ( 19951 ) we obtain a mass 
of M = 1.03 ± 0.04 Mq based on the redshift measurement 
in agreement with the mass of M = 1.10 ± 0.03 Mq based 
on the surface gravity measurement. The weighted mass av- 
erage is M — 1.075 ± 0.024 Mq, close to the ultra-massive 
range. 



3 BASED ON PARALLAX MEASUREMENTS 



Preli minary parallax information on G35-26 (ITheill et al 
Il990l) . GD 50 and PG 1658+441 was provided by iDahn 
(1999). Reading from Dahn's Figure 3, the absolute V mag- 
nitudes of these three objects are 12.9, 11.5, and 12.4, respec- 
tively. Although, the original measurements are not listed 
the absolute magnitudes may be converted into masses us- 
ing publish ed effective temperatu r es an d the mass-radius 
relations of iBenvenuto fc Althausl (ll999T l. Adopting T e s = 
43200 K for GD 50 jVennes et al.lll997b ) and T e g = 305 10K 
for PG 1658+441 (ISchmidt et all Il992n we obtain masses 
of M = 1.24 Mq and M = 1.29 Mq, respectively. These 
values are in good agreement with the spe ctroscopic masses 
l|Schmidt et al.lll992l ; IVennes et al.lll997bl ). 

Si milarly, employing t he effective temperature consid- 
ered by iTheill et alj (|l990l ). T cff = 14000 K, and the non- 



DA mass-radius relations, the absolute visual magnitude 
of G 35-26, interpolated using a grid of helium model at- 
mospheres at 7.0 ^ log<7 ^ 9.5 and logH/He = —3 and 
at a temperature of 14000 K, also suggest a high-mass of 
M = 1.16 Mq. This mass is somewhat lower th an the spec- 
troscopic mass obtained by ITheill et alj (|l990l ). M = 1.2- 
1.33 Mq. 

These c o mpari sons and the analysis of LHS 4033 



Dahnetal.1 12004) and G47-18 and ES0439-26 



(|Bergeron et all 1200 ll ) show that, although limited in 
scope, parallax measurements confirm the spectroscopic 
masses and the existence of a population of ultra-massive 
white dwarfs. 



4 BASED ON ORBITAL SOLUTIONS 

EUV surveys also uncovered a p opulation of white 
dwarfs paired with early-type stars (jVennes et alj 1 19981 ; 
iBarstow et al.l feoOll). Initial-mass to final-mass relations 
(jCatalan et al.l l2008al ) for white dwarf stars indicate that 
the progeny of early-type stars (A, B) would retain a mass 
above average and in excess of 1 Mq for progenitors with 
M > 5 Mq. The massive subluminous companion of the 
A8V star HR 8210 is such a relevant case. Far and extreme 
ultraviolet spectrosco py uncovered a DA white dwarf com- 
panion to HR 8210 (lLandsman et al.l 1 19931 ; IBarstow et al.l 
1 19941 ; I Vennes et al.|[l998l). The m ass function, f(M WD ) = 
0.219 + 0.004 |Vennes et alj|l99gj ). and the constraint placed 
on the orbital inclinati on (i < 88°) by the lack of an eclipse 
l|Landsman et al.lll995T ) restricts the mass of the white dwarf 
to M(WD) > 1.19 Mq if M(A8V) = 1.6 M . 

We confirmed the high mass for the white dwarf us- 
ing Far Ultraviolet Spectroscopic Explorer (FUSE) spectra 
that cover the white dwarf spectral energy distribution from 
« 900 A to w 1180 A (Fig. E]). The spectra (data id 
A0540909000) were obtained with FUSE on 2001 July 12 
(UT) for a total exposure time of 4199 s using the LWRS 
aperture, and were processed using CalFUS E v3.2. The 
FUSE instrumentation is briefly described by iMoos et al.l 
(2000). We fitted the Lyman line series observed with FUSE 
with sets of grid of models to constrain the effective temper- 
ature and surface gravity of the white dwarf. In total we used 
4 sets of grids, which either include or exclude the effect of 
Lyman satellite features and where line merging is treated 
either by using the formalism of Inglis fc Teller! (ll939T l (here- 
after IT) or of iHummer fc Mihalasl (Il988l) (he r eafter HH), 
which follows the treatment of iHubenv et alj ll 19941). The 
latter is described in more detail in iKawka fc Vennesl ([20061 
and lKawka et all (|2007l ). Figure [3] shows the analysis of the 
Lyman line spectra using the grid of model spectra that 
excludes Lyman satellites and used IT. The best fit model 
to the spectra corresponds to an effective temperature of 
T cff = 35490 + 70 and a surface gravity of log g = 8.85 + 0.04. 
We repeated the fit using a grid of models that include Ly- 
man satellites and use HH to obtain T a s = 34350 ± 100 
and logg = 9.04 ± 0.05. The grid of models that ex- 
clude Lyman satellites and use HH result in a best fit with 
r cff = 35000+80 and logg = 8.75+0.03. And finally the grid 
of models that include Lyman satellites and use IT result in 
a best fit with T cff = 34800 + 90 and \ogg = 9.13 + 0.04. The 
series of analyses show that the two different treatments of 
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Figure 3. Analysis of the Lyman line profiles of the DA white dwarf in HR 8210. Six separate FUSE spectroscopic channels are labelled 
and included in the analysis. Spectral ranges contaminated by geocoronal line emission were excluded from the analysis. The observed 
spectrum has been shifted by —0.5 A. 



line merging produces a difference of ~ 0.1 in the surface 
gravity, and that the inclusion of Lyman satellites increases 
the measured s urface gravity by ~ .3. We find that the best 
agreement with lVennes et all ( 199&t l is obtained using mod- 
els that include the Lyman satellites however the fit is unsat- 
isfactory as several predicted features are stro nger than ob- 
served . Similar difficulties were encountered bv lDupuis et al.l 
|2003l ) in the analysis of Lyman line profiles of the ultramas- 
sive DAp PG 1658+441. 

Em ploying the mass-radius re lations for DA white 
dwarfs (jBenvenuto fc Althausl Il999l ). we determined the 



mass and absolute magnitude of the white dwarf, M 



1.08-1.24 M© andA/v = 10.96-11.76. IVennes et atl (|l998h 
already concluded that the distance modulus derived from 
the Hipparcos parallax of the A8 star would be consistent 
with the predicted distance modulus for an ultra-massive 
white dwarf. 

To our knowledge, the case of HR8210 is rather unique 
with the possible exception of HD 2 09295. Based on their 
analysis of the 3.1 day binary period, lHandler et alj (|2002l ) 
inferred the presence of a companion to this A star. The 
measured UV excess and the constraints on the mass (M > 
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1.04 Mq) suggest that the companion is a relatively hot 
massive white dwarf, although the flux deficit in the TD- 
1 measurement at 1565 A conspires against this explana- 
tion. The case of the p utative white dwarf co mponent of 
the trip le system A Sco ( Bcrghofe r~et alj|2000h was put to 
rest by lUvtterhoeven et al.l (|2004h who inferred the pres- 
ence of a low-mass pre-main-sequence star rather than a 
hot ultra-massive white dwarf. No other ultra-massive white 
dwarfs were predicted based solely on orbital elements. It is 
worth noting that three other white dwarf companions to 
B type stars (HR 2875, 9 Hya, 1 6 Dra; IVennes et al.lll997al ; 
iBurleigh fc Barstow|[l999l. |2000| ) were uncovered from their 
extreme ultraviolet spectral signatures. Direct spectroscopic 
or orbital mass measurements are as yet unavailable. 



5 BASED ON BALMER LINE 
SPECTROSCOPY 

Fi gure 4 assem b les the spectroscopic evide nce g athered 
by iVennes et all l|l996h . IVennes et at] (|l997bh. andlVennes 



1999) based on EUV surveys and as revise d by Vennes et al. 
2008) for a total of 1 5 8 obj ects, and iKawka fc Venne" 
2006) and IKawka et all i|2007P ) based on the New Luyten 



Two- Tenths (NLTT) proper- motion catalog (66 objects). 
In the latter, using proper-motion as a proxy for distance, 
the white dwarf selection probability is primarily a func- 
tion of distance and the NLTT survey is primarily vol- 
ume limited. Due to the particular morphology of the lo- 
cal interstellar medium a nd the so-called loca l bubble (see 
iRedfield fc Linsky|[2008l '). I Vennes et all (|l997bl ) argued that 
the white dwarf selection in EUV surveys is confined to a 
volume of radius < 100 pc and characterized by low inter- 
stellar medium column density that does not hamper the 
EUV selection. They show that a white dwarf formation 
rate of 0.7 — 10 x 10~ 12 pc— 3 yr _1 is sufficient to account 
for all h ot white dwar f stars detected within 100 pc. Inci- 
dentally, ILiebert et ail (|2005h determined a formation rate 
of 0.6 x 10" 12 pc-3 yr" 1 for the DA white dwarfs in the PG 
survey. On this account, one may consider the EUV surveys 
to be volume limited as well and that few hot white dwarfs, 
irrespective of their radii (hence masses) , were missing from 
within this volume. Interestingly, the yield in ultra-massive 
white dwarfs is very similar in both NLTT, 4 out of 66, and 
EUV, 17 out of 158, surveys. A total of 21 objects out of 224 
sampled have masses in excess of 1.1 Mq, and, therefore, 
close to 10% of all white dwarfs maybe considered ultra- 
massive; 

ILiebert et ail (|2005h and IKepler et all l|2007h built 
white dwarf mass distributions based on the PG and SDSS 
surveys, respectively. Figure \5 shows the s a mple of 347 



Liebert et ail l|2005l) . Only 



DA white dwarfs analyzed by 
seven ultra-massive white dwarfs are extracted from this 
survey due to the relative faintness of these objects and 
the magnitude-limited survey strategy. Both ultra-massive 
white dwarfs hotter than 30 000 K in the PG selection are 
also part of the EUV selection. However, by applying the 
VVKnax correction to the number counts they conclude, as 
previously established in the EUV-selected count, that « 
10% of white dwarfs maybe considered ultra-m assive. A sim- 
ilar ap proach and conclusion was reached by IKepler et al.l 
J2007 j). 



o 
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Figure 4. Spectroscopically determined effective temperatures 
(T c ff) and surface gravities (logg) for the sample of EUV-selected 
DA white dwarfs (open squares) and a sample of hydrogen-rich 
DA white dwarfs from the NLTT survey (full squares; see §5). 
The measuremen ts are compared to mass-r adius relations for 
DA white dwarfs (Bcnvcnuto & Althaus 1999) for masses ranging 
from 0.3 (top) to 1.2 Mq (bottom). 
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F igure 5. Same as Fi gure [4] bur for the Palomar-Green sample 
of lLiebert et al.l l|2005l l. 



We now discuss the ultra-massive white dwarfs from the 
NLTT survey. 

5.1 WD 0457 004: NLTT 14307 

IKawka fc Vennes! (120061 ) measured a mass of 1.24 ±0.02 Mq. 
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With a temperature of 10 800 K, the white dwarf lies outside 
the ZZ Ceti instability strip but it remains useful in helping 
define the location of its red edge at high masses. 



5.2 WD 1653+256: NLTT 43827 

iKawka fc Vennesl (2006) measured a mass of 1.31±0.01 Mq. 
As in the case of NLTT 14307 and with a temperature of 
11 690 K, the white dwarf is useful in helping define the lo- 
cation of the red edge of the instability strip at high masses. 



5.3 WD 1236 495: NLTT 31372 

The star is also know as LTT 4816 and is th e most massive 
(M = 1.11 ± 0.02 M Q ; IKawka et all 120071 ) pulsating DA 
star and t he only known to dat e ultra-massive member of 
this class l|Gianninas et alj|2005l ). 

5.4 WD 1729+371: NLTT 44986 

The star is also known as GD 362 and is a peculiar 
ultra-massive (M = 1.26 ± 0.03 Mq) PAZ white dwarf 
l|Gianninas et al.ll2004l ; IKawka fc Vennes!l2006l '). The star is 
also perceived as harboring a disk of debris (Kilic et al. 
120051 ). 

However, IZuckerman et al"] l|2007l ) recently obtained 
a high-dispersion, high signal-to-noise-ratio spectrum of 
GD 362 that revealed a rich heavy-element line spectrum, 
and one or possibly two He I lines. The helium abundance 
consistent with a weak HelA5876 A line has considerable 
effect on t he effective temperature and surface gravity mea- 
surements. [Zuckerman~etj)I] (|2007T ) measured (logHe/H = 
1.1) along with a temperature T" e ff = 10540 and \ogg = 8.24. 
The mixing of hydrogen in an otherwise weakly detectable 
helium atmosphere has the effect of weakening the upper 
Balmer lines and of mimicking the effect of a high gravity. 
A parallax measurement is required to confirm the normal 
gravity for this object. 



5.5 WD 2159 754: NLTT 52728 

The star is also known as LTT 8816 and lKawka et alj (|2007i ) 
determined a mass of M = 1.17 + 0.04 Mq. This estimate is 
supporte d by the peculiar radial velocity ti(WD) = 153 ± 2 
km s _1 l|Maxted fe MarsU Il999l ) that is much larger than 
normal. To illustrate this point we let the radial velocity of 
the star vary, while fixing its proper-mot ion to the obs erved 
value, fi = 504 mas yr" 1 and 9 = 277.9° |Luvtenlll976l ). and 
compute the corresponding (U, V, W) vector. By minimizing 
the difference between this vector and the local disk vector of 
IChiba fc Beers! l|2000l ) at (U, V, W) disk = (0, -35, 0) we esti- 
mate the most probable radial velocity for the star, n ra( j = 25 
km s _1 . Therefore, the excess t)(WD) — « ra d = 128 km s _1 
represents the most likely gravitational redshift for the white 
dwarf. This redshift translates into a mass of 1.17 Mq. Of 
course, without a velocity reference point (cluster or binary 
memberships) it is not possible to confidently determine the 
gravitational redshift but it is likely to be large. 



6 SUMMARY AND CONCLUSIONS 

We critically reviewed current evidence for the existence 
of ultra-massive white dwarfs. First, we demonstrat e d that 
the high-mass white dwarfs listed in ISilvestri et alj (|200ll ) 
and presented as evidence for an ultra-massive white dwarf 
population in CPMB are in fact white dwarfs with normal 
masses (« 0.6-0.7 M ), or, in the case of G275-B16A, a 
low-metallicity F star and possible halo member. The ab- 
sence of any reliable candidates from this sample is puz- 
zling considering the large number of objects identified us- 
ing other methods. Improved radial velocity measurements 
of an enlarged sample of CPMB would help provide accurate 
gravitational redshifts and deliver the expected number of 
ultra-massive white dwarfs. 

The study of the peculiar DBAZ GD 362 exposes poten- 
tial difficulties in surface gravity, hence mass measurements 
based on hydrogen line profiles. Although spectroscopically 
evanescent at an effective temperature of 10 000 K , helium 
was found t o be the dominant consti tuent in the atmosphere 
of that star. IZuckerman et al"] (|2007T ) demonstrated that the 
reduced hydrogen abundance had the consequence of low- 
ering the surface gravity measurement to almost a normal 
level. Although helium lines are considerably stronger at 
temperatures in excess of 15 000 K, it is possible that some 
cooler ultramassive white dwarfs are in fact helium domi- 
nated but with a normal mass. 

Next, we examined the evidence based on parallax 
measurements and found that the small number of mea- 
surements, and the white dwarf masses inferred from 
these measurements are in good agreement with the spec- 
troscopic masses. Table [5] lists spectroscopically-identified 
ultra-massive white dwarfs for which parallax measurements 
are desirable. Indeed, there remain s ome questions con- 
cernin g t he validity of our app roaches, iHummer" fc Mihalad 
|l988t ) or llnglis fc Teller! (|l939h . for Balmer and Lyman line 
merging at high density and temperature. The effect of per- 
turbers on upper energy levels is essentially calibrated using 
normal gravity white dwarfs (log g = 8) and this calibration 
may not apply well at higher gravities potentially causing a 
systematic shift in mass measurements. However, a case-by- 
case review of the few spectroscopically identified high-mass 
white dwarfs with parallax measurements shows good agree- 
ment between the two methods. Table 2 lists the absolute 
visual magnitude and predicted parallax for each star based 
on parameters provided in the listed references. In addition 
to the calculated masses that we determined using the mass- 
radius relations assuming a CO core, we have also calculated 
mass estimates usi ng mass-radius relat ions for white dwarfs 
with an ONe core i|Athaus et alJ feOQS). which may be more 
appropriate for massive white dwarfs. The mass-radius rela- 
tions for ONe cores, predict masses that are systematically 
» 0.02 Mq lower than those predicted by CO mass-radius 
relations. For mas ses larger than 1.3 Mq, we u sed the mass- 
radius relations of lHamada fc Salpeterl (| 196 if ) . 

We also confirmed the high mass for the hot white dwarf 
in the binary HR 8210. A high mass was initially implied by 
the binary mass function l|Landsman et al" 1993). We fitted 
the Lyman line spectrum of the white dwarf and constrained 
the mass to be M = 1.08 — 1.24 Mq in agreement, in the 
upper mass range, with the binary parameters. 

Finally, we show that the EUV-selected population 
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of white dwarf stars is composed of ~ 10% objects with 
ma sses in excess of 1.1 Mq. A similar yield was obtained 
by iLiebert et all (120051 ) based on the PG survey and by 
iKepler et all (j2007T ) based on SDSS after large corrections 
were applied due to the magnitude-limited nature of the 
samples collected. In particular, it should be noted that only 
seven objects out of 347 from the Palomar-Green sample, or 
a fraction of 2%, met the criterion. By applying V/Vmea. cor- 
rections due to incompleteness at fainter magnitudes, the es- 
timated fraction was re-evaluated at 10% in agreement with 
the yield directly measured in the EUV selection. 

The origin of ultra-massive white dwar fs remains un- 
certain . Initial-mass to final-mass relations l|Catalan et al.l 
2008a) indicate that main sequence stars with masses in 
excess of « 6 Mq generate white dwarfs with masses in ex- 
cess 1 Mq , a situation best illustrated by the massive white 
dwarf (LB 1497) me mber of the Pl e aides. By re-evaluating 
available cluster data lCatalan et all (|2008al ) revised the final 
masses upward, and managed to reproduce the high-mass 
peak in both SDSS and PG empirical mass distributions. It 
is therefore possible that white dwarfs with masses in excess 
of 1.1 Mq are the products of single star evolution and that 
the binary merger scenario may only apply to a minority 
of peculiar objects such a s the fast rotating magnetic white 
dwarf WD 0325-857 (see I Vennes et al.ll2003l . and references 
therein). 

The existence of a substantial population of ultra- 
massive white dwarfs supports the concept of a steeper 
initial-mass to final-mass relations linking 6 Mq progenitors 
with > 1.1 Mq white dwarfs as proposed bv lCatalan et all 
l|2008al ). Ultra- massive white dwarfs in close binaries are also 
likely Type la supernova progenitors (jParthasarathv et al.l 
120071 ). 
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Table 2. NLTT and EUV spectroscopic sample of ultra-massive white dwarfs 



WD 


Other names 


Mass" 


M v 


Mass 6 


Predicted it 


Ref. 








(M@) 


(mag) 


(M & ) 


(mas) 




0457- 


-004 


NLTT 14307 


1.24 ±0.02 


14.09 


1.22 ± 0.02 


58.8 


1 


1236- 


-495 


NLTT 31372 


1.11 ± 0.02 


13.13 


1.09 ±0.02 


66.7 


2 


1653- 


f256 


NLTT 43827 


1.31 ± 0.01 


14.29 


1.28 ±0.01 


30.3 


1 


17294 


-371 c 


NLTT 44986 


1.26 ±0.03 


14.54 


1.24 ±0.03 


45.5 


1 


2159- 


-754 


NLTT 52728 


i i 7 +0.06 
1 - 1 '-0.07 


14.30 


1 14+0.05 
1 - 1 ^-0.06 


71.4 


2 


0001- 


f433 


EUVEJ0003+435 


1.18 ±0.04 


11.19 


1 16+ 004 
1±D -0.05 


7.2 


3 


0136+251 


EUVEJ0138+-253, PG 


i io+0.04 
1 " 1O -0.05 


10.97 


1.10 ±0.05 


9.9 


3 


0325- 


-857 


EUVEJ0317-855 


1.32 ± 0.03 


12.64 d 


1.32 ±0.03 


37.0 d 


2, 4 


0340- 


f-103 


1RXSJ0343+1030 


i 99 +0.04 
X - -0.05 


11.01 


1 20+ 05 

i ' ZU -0.06 


8.0 


3 


0346- 


-Oil 


EUVEJ0348-009, GD 50 


1.23 ± 0.02 


11.49 


1.21 ±0.02 


31.3 


3 


0440- 


-038 


EUVEJ0443-037 


1.14 ±0.07 


10.57 


1.11 ±0.07 


5.4 


3 


0554- 


-165 


1RXSJ0557-1635 


1 ori+0.07 
1 "5U_ 12 


11.86 


1 oq+0.09 
i ' 2S -0.12 


5.4 


3 


0630- 


f200 


1RXSJ0633+2001 


i 9O+0.04 
1,zz -0.05 


11.11 


i on 4 " 0,05 


6.9 


3 


0652- 


-563 


T^T TVFT Tnf^'i ^fizl 


i lfi +0.05 




i iQin 
i . ±o m u.ud 


inn 


Q 



0821- 


-252 


1RXSJ0823-2525 


1.27 ±0.05 


11.83 


i 9C-+0.04 
1 - zo -0.05 


12.2 


3 


0914- 


-195 


EUVEJ0916-197 


-, 9 7 +0.06 
± " a -0.07 


11.47 


i 9C+O.O6 
1 -^°-0.08 


6.8 


3 


1022- 


-301 


1RXSJ1024-3021 


1.21 ± 0.04 


11.63 


1.19 ±0.05 


13.3 


3 


1529- 


-772 


EUVEJ1535-774 


1 26+ 03 
1 ' ZD -0.04 


11.54 


i 94+0.03 
-0.04 


10.6 


3 


1658- 


(-441 


EUVEJ1659+440, PG 


1 34+O.O6 
1 -° -0.14 


12.83 


-, oq+0.07 
±.oo_ Q 16 


43.5 


3 


1711- 


f668 


1RXSJ1711+6645 


i 99+O.O3 
X - -0.05 


11.25 


1 20+ ' 04 
i - ZU -0.05 


6.8 


3 


1724- 


-359 


EUVEJ1727-360 


1 20+ - 02 
±zu -0.04 


11.66 


1.17 ±0.03 


17.2 


3 


1740- 


-706 


EUVEJ1746-706 


1.17 ±0.04 


11.01 


1.14 ±0.04 


7.9 


3 



a Using the mass-radius relations for a CO core (Bcnvcnuto & Althaus 1999). For 
masses above 1.2Mq, mass-radius relations of lHama da &: Salpeterl lll96lfl were used. 

Using the mass-radius relations for a ONe core llAthaus et al.|l2005h . For masses 
above 1.3Mq, mass-radius relations of lHamada &; Salpete rl lll96lh were used. 
c Based on the model fit of Kawka & Vennes (2006) ass uming He/H = 0. d Based 

on the distance modulus of the DA com panio n LB 9802 iKawha et aLll2007f). 

Refere nces: (!) iKawka fc Vennesl ||200i); (2) iKawka et al.l i2007f) ; (3) IVennes et alj 
d2008t) : (41 1 Vennes et al.1 ll2003h 
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